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MONCOR — A  Program 
to  Compute  Concordant  and  Other 
Monotone  Correlations 


by 


George  Kimeldorf,  Jerrold  H.  May,  and  Allan  R.  Sampson 
University  of  Texas  at  Dallas,  University  of  Pittsburgh,  and 
University  of  Pittsburgh 


ABSTRACT 

’ K  User's  Manual  is  provided  for  the  interactive  FORTRAN  program  MONCOR 
Which  coalites  the  concordant  monotone  correlation,  discordant  monotone  corre¬ 
lation,  isoconcordant  monotone  correlation,  isodiscordant  monotone  correlation 
and  their  associated  monotone  variables.  Data  input  are  ordinal  contingency 
tables  in  matrix  form^ _ 
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USER'S  MANUAL:  MONCOR-- 


A  Program  to  Compute 

Concordant  and  Other  Monotone  Correlations 

by 

1  2 
George  Kimeldorf  ,  Jerrold  H.  May,  and  Allan  R.  Sampson 

University  of  Texas  at  Dallas,  University  of  Pittsburgh,  and 

University  of  Pittsburgh 

V 

1.  General  Background 

MGNCOR  is  a  set  of  routines  to  compute  monotone  correlation  measures  for 
two  dimensional  probability  distributions  of  random  variables  X  and  Y  given  by 
matrices  of  the  form 


Pll 

P12 

’•*  P1J 

P21 

• 

• 

P22 

• 

• 

’**•  P2J 

• 

Pll 

• 

PI2 

...  pTJ 

where  p^.  *  ProbfX  *  x^,  Y  ■  y^),  i  ■  1,  ...,  I  and  j  =  1,  ...,  J.  These  proba¬ 
bilities  may  be  theoretical  probability  distributions  or  may  be  estimates  derived 
from  ordinal  contingency  tables. 

The  general  approach  followed  in  MONCOR  is  based  upon  the  algorithm  given 
in  Kimeldorf,  May  and  Sampson  (1980)  with  a  more  theoretical  development  in 

*The  work  of  this  author  was  supposed  by  the  National  Science  Foundation 
under  Grant  MCS-8002152. 
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Kimeldorf  and  Sampson  (1978) .  MONCOR  computes  the  concordant  and  discordant  mono¬ 
tone  correlations,  and  their  isoscale  versions,  by  formulating  a  linearly  con¬ 
strained  nonlinear  programming  problem  whose  optimum  value  is  the  desired  monotone 
correlation  measure.  A  specially  modified  form  of  the  nonderivative  nonlinear 
programming  algorithm  QRMNEW  (see  May  (1979))  is  used  to  solve  the  nonlinear 
program. 

2.  Program  Structure  and  Coding  Considerations 

2.1  Structure 

MONCOR  consists  of  two  basic  components:  routines  specific  to  monotone 
correlation  and  routines  adapted  from  the  QRMNEW  program.  A  brief  overview  of 
the  calling  sequences  is  given  in  Figure  2.1.  The  primary  goal  throughout  the 
coding  of  MONCOR  was  simplicity  of  presentation.  The  components  of  QRMNEW 
adapted  for  MONCOR  are  written  in  part  in  a  non -structured  style,  but  with  read¬ 
ability  enhanced  by  comment  lines. 

Since  MONCOR  was  designed  to  be  used  in  an  interactive  mode  for  small  prob¬ 
lems,  dense  matrix  storage  is  utilized  throughout.  Array  storage  is  made  available 
to  all  subroutines  through  labelled  common  only.  All  dimensions  are  explicitly 
set  in  each  subroutine.  MONCOR  version  2.3  is  dimensioned  so  as  to  handle  proba¬ 
bility  matrices  of  dimension  up  to  20  x  20.  The  source  listing  comments  indicate 
how  to  redimension  the  arrays  for  larger  problems. 

The  original  source  code  for  MONCOR  consists  of  1932  statements  written  in  DEC 
FORTRAN-10  and  is  implemented  on  the  University  of  Pittsburgh's  DEC  1099. 

2.2  Programming  Considerations 

(a)  Types.  Double  precision  arithmetic  is  used  throughout,  and  both  variables 
and  arrays  are  set  by  an  IMPLICIT  statement  at  the  beginning  of  each  routine. 
Without  exception,  the  type  of  a  variable  is  given  by  the  first  charactei  of  its 
name,  as  follows: 


SfiAROIf 


Calling  Structure  For  MONCOR  Routines 


A-H,  0-Z 


Double  Precision 


I-N  Integer. 

No  half  word  integer  (INTEGER  *2)  variables  are  used.  No  logical  variables  appear 
in  MONCOR. 

(b)  COMMON.  Blank  COMMON  is  not  used.  There  are  21  labelled  common  blocks, 
containing  all  the  arrays  and  important  scalar  variables. 

(c)  EQUIVALENCE.  No  EQUIVALENCE  statements  are  used. 

t 

(d)  DATA.  No  DATA  statements  are  used. 

(e)  Arrays-.  There  are  both  one-dimensional  and  two-dimensional  arrays,  none  of 
variable  length.  Declarations  of  the  form  REAL  A(l)  are  never  used. 

(f)  Subscripts.  Subscript  expressions  of  the  form  X(M  +  J)  as  well  as  X (INDEX  (J)) 
are  often  used. 

(g)  Avoidance  of  Trouble  Spots.  The  following  FORTRAN  syntactical  items  are 
not  used: 

Arithmetic  IF;  Assigned  GO  TO;  BACKSPACE;  Blank  COMMON;  BLOCK  DATA; 
COMPLEX;  ENDFILE;  ENTRY;  EQUIVALENCE;  EXTERNAL;  FIND;  FUNCTION;  NAME- 
LIST;  Nonstandard  RETURN;  PAUSE;  PRINT;  PUNCH;  TIME;  Uncommitted  READ 

and  WRITE. 

2.3  Conversion  of  MONCOR  to  Other  Machines 
If  double  precision  is  not  required,  the  IMPLICIT  statements  should  be 
deleted,  and  the  names  of  the  library  functions  (DABS,  DMAXl,  DMIN1 ,  DLOG,  DSQRT, 
etc.)  should  be  changed  to  their  real  equivalents  (ABS,  AMAX1 ,  AMIN1 ,  ALOG,  SQRT, 
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If  double  precision  is  required  but  the  IMPLICIT  statement  is  not  allowed, 
all  variables  starting  with  A-H  and  0-Z  must  be  typed  explicitly.  Such  a  con¬ 
version  will  require  many  lines  of  otherwise  unnecessary  code,  and  will  also  be 
prone  to  error,  since  FORTRAN  compilers  do  not  usually  issue  warnings  for 
undeclared  variables. 

In  FORMAT  statements,  apostrophes  are  used  to  delimit  strings.  Certain 
strings  in  the  main  program  contain  the  character  *. 

The  main  program  calls  the  FORTRAN-10  basic  external  function  RAN,  which 
returns  a  [0,1]  random  number.  The  call  to  RAN  may  be  modified  to  accomodate  any 
other  random  number  generator. 

3.  Program  Description 

In  all  sample  M0NCOR  outputs,  those  items  inputted  by  the  user  are  underlined. 
On  this  output,  references  to  material  in  Section  3  are  typed  in  the  left  margin. 

3.1  Example 

To  describe  the  use  of  MONCOR,  the  British  Mobility  data  example  is  employed. 
(See  Kimeldorf,  May  and  Sampson;  the  original  data  are  in  Glass  and  Hall  (19S4, 
p.  183).)  The  estimated  probability  matrix  (p^j)  derived  from  Table  5.1  of  Kimel¬ 
dorf,  May  and  Sampson  is  given  in  Table  3.1. 

Table  3. 1 


Estimated  Probability 

Matrix 

British  Mobility  Data 

Son's  Occupational 

Status 

S2 

S3 

S4 

S5 

Si  .014 

.013 

.002 

.005 

.002 

Father's  Occupational 

S2  .008 

.050 

.024 

.044 

.016 

Status 

S3  .003 

.022 

.031 

.064 

.027 

S4  .004 

.043 

.053 

.205 

.128 

S5  .001 

.012 

.021 

.091 

.117 

In  this  example  I  »  J  « 

5. 

5 


3, 2  Data  Input 

The  program  requires  input  in  the  following  order: 

(a)  the  number  of  rows  of  Ip.^}: 

(b)  the  number  of  columns  of  (p. .}:  J, 

i  j 

(c)  the  device  number  that  the  matrix  is  to  be  read  in  from  (entering 


a  5  allows  reading  from  terminal). 


(d)  the  values  for  {p*j)  read  in  a  row  at  a  time  (it  is  assumed 
max(I, J)  <  20). 

After  {Pjj)  is  inputted, 

(e)  the  matrix  is  printed  by  M0NC0R  under  the  heading  INPUT  MATRIX. 


An  example  of  this  procedure  is  given  for  the  British  Mobility  data. 

3.2(a)  HOW  MANY  ROUS  IN  YOUR  MATRIX?  5 

3.2(b)  HOU  MANY  CfHUMNS  IN  YOUR  MAIKIX?  ~ 

3.2(c)  ON  WHAI  DtVICr  IB  YOUR  MATRIX?  I.N1FR  f~ IT  YOU  WISH  10  IYII  II  IN  Nil W.  !, 

ENTER  THE  MATRIX  A  ROW  A!  A  I  I  ME-  MAX  20  NOS.  Ilk  I  INI 


>  ,014 

.013 

,002 

.005 

.002 

>.000 

.  olio 

#  644 

TtflTT 

,003 

.022 

.031 

•  0A4 

.027 

>  ,  004__ 

.043 

.053 

.205 

.  1 2H 

;  ]ooi 

f o  1 2 

.021 

.<)<?! 

.  it? 

rwptiT 

MATffl  X 

.0140 

,0130 

.00 

20 

.  0050 

,  0020 

.0080 

.  0500 

.02 

40 

.0440 

.01  AO 

2fc)  .0030 

.  0220 

.03 

10 

•  0A4»> 

.  02  70 

.0040 

.0430 

.  05 

30 

.  2050 

,  1  2U0 

*0010 

,0120 

.02 

10 

.0910 

.1170 

r  Nn  m 

IMH.II 

3 . 3  Output 

(a)  After  the  matrix  (p.j)  has  been  output  by  MONCOR,  the  user  must  specify 
which  one  measure  of  monotone  correlation  is  desired.  The  user  needs  only  to 
enter  one  number  (1,  2,  3  or  4)  according  to  the  MONCOR  prompts. 

(b)  For  each  monotone  measure  calculation,  the  user  has  the  option  of 

calling  for 

(i)  intermediate  point  output  during  the  optimization  process 
(ii)  inputting  a  starting  point  for  the  optimization  process. 


Responses  to  MONCOR  questions  concerning  (i)  or  (ii)  arc  YP.S  or  NO. 
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For  the  remainder  of  the  output  description  in  Subsection  3.3,  it  is  assumed  that 
neither  option  (i)  nor  (ii)  is  requested. 

(c)  MONCOR  then  evaluates  the  correlation  at  the  monotone  extreme  points 
(see  Section  4,  Kimeldorf,  May  and  Sampson).  The  number  of  such  monotone  extreme 
point;  u  v  n  by  MONCOR  (if  the  concordant  or  discordant  monotone  correlation 
is  requested,  the  number  of  monotone  extreme  points  is  (I-1)(J-1);  if  the  iso- 
concordant  or  isodiscordant  monotone  correlation  is  requested,  the  number  of 
monotone  extreme  points  is  (1-1)). 

* 

(d)  MONCOR  computes  the  largest  correlation  among  these  monotone  extreme 

points,  and 

(e)  prints  one  monotone  extreme  point  at  which  this  maximum  is  achieved. 

(f)  MONCOR  then  prompts  the  user  to  enter  a  five  digit  number  which  is  used 
to  warm  up  the  FORTRAN-10  basic  external  function  random  number  generator  RAN.  RAN 
is  called  that  number  of  times  before  being  used  to  generate  the  ten  random 
starting  points  (see  Kimeldorf,  May  and  Sampson  (Section  41). 

(g)  The  estimated  monotone  correlation  measure  is  given  based  upon  the 
10  random  starting  point  calculations  as  well  as  the  monotone  extreme  points 
(Kimeldorf,  May  and  Sampson  (Section  4)). 

(h)  A  point  (monotone  variable)  at  which  the  monotone  correlation  is 
achieved  is  printed  by  MONCOR.  MONCOR,  without  loss  of  generality,  sets  X(l)  = 

Y(l)  =  0  and  X(I)  -  Y(J)  *  1.  If  an  iso-measure  is  selected,  only  the  X  mono¬ 
tone  variable  is  printed.  (The  Y  monotone  variable  is  identical  by  definition.) 

(i)  MONCOR  then  asks  if  the  user  wants  to  use  the  inputted  {p. matrix 
to  calculate  other  measures  of  monotone  correlation.  If  the  user  responds  YF.S , 
it  begins  at  3.3(a)  again.  If  the  user  responds  NO,  an  opportunity  is  given  to 
input  another  {pjj}*  (see  3.2(a))  or  to  terminate  MONCOR. 

For  illustration,  the  calculation  of  the  isoconcordant  monotone  correlation 
for  the  British  Mobility  data  is  given  below. 


7 


5.3(a)  MIHUH  ONE  OF  IHI  t  Ml  I  OWING  DO  VIHI  WISH  III  Sfl  * 

i  riiNcnkUANi  hiihiiiiini  i.orkf i  a i  ion 
?-nt!iU»RriANi  HiiMiiriim  cim-rf  i  amon 
j  innroNriimANi  mono i  oni  uiruiai  ion 
4- isoinsriiKiiANi  wimmim  cokri  i  ai  ion 


**************** tt ** **+tt  ****  *********  ******  t  ‘UK  MINI  'ORMAN  I  MONO  I  ONI  (  I  ll-.'Kt  I  A  I  Ini)  *♦*+♦! 
3.3(b)(i)  no  YOU  WISH  HI  Off  I  NURMI  MAIF  point  atirrui?  NU 

3.  3(H)  ( i  i)  I'M  YOU  WISH  10  I  Will  A  GIARIING  POINT*1  NO 

3.3(c)  f«ASfn  UN  4  MONO  ION)  E  X  IRFMF.  POINTS.  I  ML 

3.3(d)  MAX  I  Ml  IN  Or  P'OKRFI  A I  KINS  AMONG  MONOTONE  CXTRIMF  POINTS  IS  <> .  P>.  .V'.  r -(»  0  1 

3.3(c)  X<  1 )  -~  .00000  X<  2>- 1.00000  X<  3>~1.  00000  X<  4>  1.00000  >.(  *,)  -  I  .0«'i>nn 


3.3(f)  rNTER  A  fi  IUG1I  INIFGFR  40147 

3.3(r)  Ff*  1  IMAIRli  rsnCONrORUANT  MONO  ION!  rORRfll  A)  ION  IS  0 . 4<P'.0.HV,  7  Al 

3.3(h)  X'<  l)^  ,00000  X(  :*»>  .A?7/H  X(  3>-  ,  17404  4  y<  4>  .  W4(*V  y<  %  ^ 

3.3(i)  ANOTHER  RUN  Utlll  IMF  SAMI  MA1K1X?  VI  ^ 


3,4  Option:  Inputting  a  Starting  Point 

Unless  otherwise  specified,  M0NCOR  generates  10  random  starting  points  from 


which  the  optimization  procedure  starts.  However,  the  user  has  the  option  of 
inputting  a  monotone  starting  point  by  replying  YES  to  the  question  "PO  YOU  WISH  TO 
INPUT  A  STARTING  POINT?"  The  user  must  then  input  X(2) ,  X(T-l)  and  Y(2) , 

Y(J-l).  (If  an  iso-measure  is  requested,  then  only  X(2) ,  .  ..,  X(I-l)  must  be 

input.)  MONCOR  sets  S(l)  =  Y(l)  =  0  and  X(I)  =  Y(J)  =1.  If  this  option  is 
employed,  no  monotone  extreme  points  are  evaluated  (see  3.3(c))  nncl  no  random 
starting  points  are  used  (see  3.3(g),  3.3(f)). 

For  the  British  Mobility  data,  the  following  illustrates  inputting  a  pair 


of  starting  points  for  the  concordant  monotone  correlation. 


00  YOU  WISH  10  INPUT  A  STARTING  PUtMT* 
X< 

x<  3) ‘ 

3.4  X<  4)— — 

Y  (  ? ) 

Y  <  3) 

Y  (  4  ) ;-’ 


3.5  Option:  Requesting  Intermediate  Point  Output 

Unless  otherwise  specified  MONCOR  does  not  print  out  the  intermediate  steps 
in  the  iterations  of  the  nonlinear  optimization  (see  Section  4  of  Kimeldorf, 

May  and  Sampson).  However,  the  user  has  the  option  of  requesting  intermediate 
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output  by  replying  YES  to  the  question  MD0  YOU  WISH  TO  SEE  INTERMEDIATE  POINT 
OUTPUT?” 

If  the  user  has  not  inputted  a  starting  point,  MONCOR  will  print  the  inter¬ 
mediate  evaluations  for  each  of  the  10  random  starting  points.  If  the  user  has 
inputted  a  starting  point,  MONCOR  will  print  out  the  intermediate  evaluations 
based  upon  the  given  starting  point. 

For  each  iteration,  MONCOR  outputs 

(a)  the  iteration  number  (with  iteration k number  1  corresponding  to  the 
random  or  inputted  starting  point) , 

(b)  the  number  of  evaluations,  which  is  the  number  of  times  QRMNEW  has 
evaluated  the  objective  function  cumulatively  up  to  the  given  itera¬ 
tion. 

(c)  the  correlation  based  upon  the  current  iteration, 

(d)  the  gradient  which  measures  the  degree  of  the  local  steepness  at 
that  point  (see  a  description  of  QRMNEW  in  May  for  details), 

(e)  the  coordinates  of  the  current  point  (X(l) ,  X(I),  Y(l),  Y(J)  for  the 

current  point  are  suppressed,  and  for  iso-measures,  Y(2),  Y(J-l) 

are  not  printed) . 

(f)  When  the  algorithm  stops,  MONCOR  prints  "OPTIMAL  SOLUTION  FOUND"  and 
provides  the  optimizing  point,  the  value  for  the  final  iteration, 
and  the  total  cumulative  number  of  function  evaluations  (see  3.5(b)). 

Again  for  the  British  Mobility  Data,  the  intermediate  point  output  for  the 
isoconcordant  monotone  correlation  is  given  below  for  random  point  number  2  (no 
starting  point  was  input) . 


RANDOM  rniNT  MtlMIM  f 


3.S(i) 

ITN 

FVAL 

CORR 

ORAP 

X(  7) 

X  <  3  ) 

X(  4) 

3.5(h)  1 

1 

10 

0.4597 

3. HP  01 

.  529 

.444 

.  7  75 

3.5(c)  j 

n 

?0 

0.4794 

1 . 7P  Ot 

.7  74 

.9  77 

1  .000 

3.5(d)  j 

3 

30 

0. 4 96 A 

3 .  on  o;» 

.  a;»a 

.034 

.913 

3.50)  | 

4 

40 

0 » 4940 

7.  ad  04 

.a;*b 

.044 

.974 

,  xxxxxxxx 

nn  FMAL  SOI  IK  ION  FOUND 

3.5(f)  j 

x« 

1  >*■  , 

,00000 

xc  :»> 

.42770 

X(  3 

o-  .n 

t  .O-'Of'O  V* 


CORRFl  ATKIN  OOF FF IC I TN I - 


0.49AH21 


I  I  INI  J  1 1  IN  I  VAI  1 1 A  I  IIJM*; 
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3.6  Cogent 

A  User  who  wishes  to  compute  the  (Pearson)  correlation  for  a  given  set  of 
values  of  X(2),  X(I-l)  and  Y(2),  ...»  Y(J-l)  can  use  Options  3.4  and  3.5, 

so  that  the  first  iteration  evaluation  will  provide  the  desired  correlation. 

For  example,  to  compute  Spearman's  rank  correlation,  set  X(i)  =  (i-l)/(I-l)  for 
i  =  2,  ....  1-1,  and  Y(j)  =  (j-l)/(J-l)  for  j  =  2,  ....  J-l. 
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